Zero Droplet Stiffness Exponent is Revealed in Short Range Spin Glasses when 
Probed with Large Avalanches Induced by Long Range Interactions 
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We probe the droplet excitations in short range spin glasses by adding a perturbative long range 
interaction that decays with distance as a power law: Jjr" . It is shown that if the power law 
exponent cr is smaller than the spatial dimension d, the perturbation induces large scale avalanches 
which roll until they force the system to develop a pseudo gap in the excitation spectrum of the 
stabilities. This makes the perturbative long range interactions relevant for a < Oc = d. The droplet 
theory predicts that the critical exponent ctc depends on the droplet stiffness exponent a& = d — Q . 
Combining these two results leads to a zero stiffness exponent ^ = in the droplet theory of short 
range spin glasses. 
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Spin glasses, and frustrated models in general, are rel- 
evant to a mesmerizingly wide variety of problems, such 
as random magnets, disordered electronie systems, neural 
networks, and even biological pathways. In spite of the 
pressing need for an understanding of their physics, even 
after more than thirty years of study, the community re- 
mains profoundly split regarding the correct theoretical 
framework for short range spin glasses. The two main 
approaches are the replica symmetry breaking (RSB) 
%, S, H H and the droplet pictures. [1, d, 0, Q These 
theories can be distinguished, for example, by their pre- 
dictions for the energy E[V) required to reverse a cluster 
or droplet of spins characterized by a linear size L. In 
the spin glass phase, the scaling form E[V) ~ defines 
the droplets stiffness exponent d, which is positive in the 
droplet theory: Q > but zero in the RSB theory: 6' = 0. 

Whether the stiffness exponent is finite or not, has pro- 
found consequences for the essential physical mechanism 
of the spin glass transition. In the droplet picture the 
driving force of the glass transition is the emergence of 
the finite stiffness {Q > 0) of the droplets. In contrast, in 
the RSB picture the stiffness remains weak {9 = 0) across 
the glass transition, which is driven instead by another 
mechanism: the divergence of barriers which fragment 
the free energy landscape. 

There are many other differences between the two the- 
ories predictions, such as whether the droplets are fractal 
or not. A recent attempt to bridge those differences in- 
troduced a hybrid approach, called the trivial-non-trivial 
(TNT) theory [i, [la, [HI, El, E, El , which is built on the 
notion that different types of excitations of the same size 
L may have different associated energy costs. Droplet 
excitations with size L cost an energy Edropiet{L) ^ 
whereas domain wall excitations of similar size cost an- 
other energy Edomain —waiiiL) ^ L"''™ . Direct numeri- 
cal measurements by the creators of the TNT approach 
demonstrated that the two stiffness exponents are differ- 
ent for the short-ranged Edwards- Anderson spin glass in 
3d: the droplet stiffness exponent 9 is indistinguishable 



from zero, whereas the domain wall exponent is clearly 

positive flrf^ > 0. 0, E3, EH, El, El, El 

However, all numerical works could be criticized for 
having been performed on systems with unconvincingly 
small sizes. Thus, numerics is sometimes assigned sec- 
ondary importance in the quest for the ultimate spin glass 
theory. [H 

The present Letter aims to contribute to this discus- 
sion by adopting the theoretical framework of the droplet 
theory and deriving a theoretical bound for the droplet 
stiffness exponent 9 purely on theoretical grounds with- 
out resorting to numerics. We probe the ground state 
of a short range spin glass by adding a long range inter- 
action with a small amplitude, investigate the relevance 
of this perturbation, and compare our result with the 
considerations of the original droplet paper. [^ In par- 
ticular, we determine how the distribution of local sta- 
bilities P{X) changes when the perturbative long range 
interaction is introduced. We find that the perturbation 
induces avalanches which open a pseudo gap in the sta- 
bility distribution P(A) and therefore become relevant if 
the exponent of the power law interaction a is below its 
critical value CTc, equal to the spatial dimension d. 

The original droplet paper also determined the critical 
value of the a exponent and found that CTc depends not 
only on the spatial dimension d, but also on the stiffness 
exponent of the droplets: ac = d — 9. The combina- 
tion of these two exponent relations yield a zero droplet 
stiffness exponent, 9 — 0, purely on analytical grounds, 
which is our central result. 

In spin glasses, the local dynamics and energetics can 
be conveniently characterized in terms of stabilities which 
are the products of the local spins Si and their respective 
local fields Hi : Xi = SiHi = Si J2j Jij^j- At zero tem- 
perature, all spins have positive stabilities in any locally 
stable state, including the ground state, and the values 
of these local stabilities describe the energy of one-spin 
excitations from that stable state. For short range inter- 
actions, the density P{X) of these stabilities assumes a 



finite value at zero stability: P(A = 0) > 0.[l3| Thus, a 
fraction of the spins experience vanishing effective fields 
even in the ground state. As such, they may be expected 
to be particularly sensitive to perturbations. 

The energetics of the short range system is tested by 
adding perturbations characterized by weak random long 
range interactions whose amplitude decays with distance 
as a power law with an exponent tr: J/ r°^- , J being a ran- 
dom number with a Gaussian distribution of mean and 
width J.[l3| We evaluate the relevance of this perturba- 
tion at zero temperature by analyzing whether the spin 
glass ground state is stable against avalanches induced 
by this perturbation. 

The stability analysis starts by flipping a single spin 
after having added the long range perturbation and de- 
termining the average number of spins Nonstable which 
get destabilized by this single flip. Under a sequential 
spin-updating dynamics, one of these new unstable spins 
will flip in the next time step. Thus, if a single spin flip 
destabilizes more than one additional spin on average, 
i.e., Nunstabie > 1, then the spin flip starts an avalanche, 
or exploding chain reaction, causing a large scale desta- 
bilization of the ground state of the unperturbed short 
range spin glass. If less than one spin is destabilized, i.e., 
Nunstabie < 1, the avalauchc dies out quickly, indicat- 
ing that the ground state was stable against the pertur- 
bation. Therefore, the avalanche stability criterion lies 
at exactly one newly destabilized spin per spin flip, i.e., 

^unstable — f- 

Fig. 1 illustrates the calculation of Nunstabie- Flipping 
a spin at site i can destabilize spins at a distance r from 
the site i if the stabilities of the latter spins A(r) arc 
smaller than their interaction with the flipped spin at 
site i: A(r) < J{r) = J/r°' . For a short range spin glass, 
the stability distribution P{X) can be approximated at 
small stabilities by its value at zero stability -P(O) > 0, 
yielding the following estimate for the average number of 
newly destabilized spins: 

Nunstable{r) = / d\P{\) « J(r)P(0) = P(0)J/r'^ 



NunstableXR) ^ drr^-' Nunstable{r) ^ R''-" (2) 

Jo 

Here, Eq. (2) shows the summation of the unstable spins 
for distances up to the system size R. This estimate dis- 
plays our central result clearly: for a < d, a flip of a sin- 
gle spin destabilizes a large number of spins via the long 
range perturbative interaction: N^nstabieiR, o' < d) ~ 
jld-cr >> 1^ irrespective of how small the amplitude of 
the perturbation. This is so, because even though the 
strength of the long range interaction J{r) diminishes 
with increasing distance r, but this effect is more than 
compensated by the number of spins impacted by J(r) 
increasing with r, if cr < d. The subsequent flipping of 
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FIG. 1: Tlie calculation of N^nstabieir) ■ Upon the flipping 
of a single spin, tlie spins at distance r witfi stabilities A < 
J(r) become unstable. Here P(A) is the distribution of the 
stabilities A for the short range spin glass 

these destabilized spins destabilizes an even more popu- 
lous generation of spins, setting off a large avalanche. 

At the start of the avalanche, is it appropriate to an- 
alyze the number of spins, destabilized in subsequent 
steps, independently since the spatial structure of long 
range perturbation is not correlated in any way with the 
ground state configuration, created by the short range in- 
teraction. While the rolling of the avalanches may start 
to develop such correlations, these should be of secondary 
importance as each spin flip destabilizes a large number 
of spins, thus sustaining the avalanche. 

Similar to the Efros-Shklovskii scheme for Coulomb 
glasses [3], one expects that the avalanches will keep 
rolling until the system qualitatively transforms itself and 
eliminates the fuse of the avalanche, which is the large 
number of spins with low stabilities. Formally this is 
represented by P(0) being greater than for the short 
range spin glass. The avalanches will roll until the av- 
erage number of the destabilized spins Nunstabie is sup- 
pressed to or below 1. The minimal transformation of 
the system by which this stabilization can be achieved is 
that P{X) develops a power law form at small A stabil- 
ities: P(A) ~ A". In the Coulomb Glass literature this 
process is called the opening a pseudo gap or soft gap. 
Indeed, a positive exponent a of the stability distribution 
reduces the number of destabilized spins: 

Nunstabieir) = / rfAP(A) « J(r)"+i (3) 
Jo 

NunstableiR) = / drr''-^ Nunstable{r) ~ P'*--(" + l) 

Jo 

(4) 

From Eq.(4), this power law Ansatz for P(A) is capable 
of arresting the avalanches and stabilizing the system for 
a > ac — {d/a) — 1 by making the number of destabi- 
lized spins Nunstabie{R) non-divcrgcnt with the system 
size R. The corresponding opening of the pseudo gap 
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FIG. 2: The P(A) distribution of stabilities develops a P(A) ~ 
A" soft gap driven by the avalanches induced by the pertur- 
bative long range interaction. 



in P(A) is illustrated in Fig. 2. Following the insights 
from the Coulomb gap literature, one expects that the 
system develops these large scale rearrangements only to 
the minimal degree necessary for stopping the avalanches. 
In other words, the system drives itself only to the state 
of marginal, or critical stability. This is achieved by the 
exponent a assuming the critical value ac- 



a = etc 



[d/a) - 1 



(5) 



which saturates the Efros-Shklovskii bound. In this 
sense, systems with long range interactions are similar 
to those which exhibit Self Organized Criticality (SOC). 
Note that Eq.(5) holds only in the regime d/2 < a < d. 
For (T < d/2 the interaction has to be rescaled to en- 
sure that the energy remains an extensive quantity. With 
the proper rescaling of the interaction a remains pegged 
at a{d/2) ~ 1 for < d/2, including the Sherrington- 
Kirkpatrick case of cr = 0. For a > d the system re- 
mains effectively short ranged. Accordingly, no large 
scale avalanches form and no power law soft gap is needed 
for stability. 

In our earlier work on the hysteresis in the Sherrington- 
Kirkpatrick model we observed and reported the forma- 
tion of such a pseudo gap with an exponent of rn^st = 1 
for the states reached by a hysteretic history, [l^ This 
exponent ahyst = 1 is entirely analogous to the above 
a{a = 0) = 1 result for the ground state. In the work 
on hysteresis we also interpreted the saturation of the 
exponent bound as a sign of the system driving itself es- 
sentially to a Self Organized Critical state. 

Returning to the central inquiry of the present work, 
Eq.(5) expresses the same criterion as Eq.(2) in the fol- 
lowing sense. For a < ac = d, a becomes greater than 
zero and correspondingly the system is forced to respond 
to the long range perturbation by large scale avalanches, 
qualitatively rearranging the ground state and develop- 
ing a pseudogap in P(A). In other words, for cr < CTc = d 
the long range interaction is a relevant perturbation. 

This relevant perturbation should be detectable 



through, for example, the low temperature behavior of 
many measurable quantities, such as the specific heat, 
the magnetic susceptibility, or, in electronic equivalents, 
the tunneling conductance. 

To complete our argument, we presented two separate 
analysis of the impact of a long range perturbation on a 
short range spin glass, based on (i) computing the num- 
ber of unstable spins in avalanches induced by the per- 
turbation, and on (ii) describing the development of a 
pseudogap, induced by the perturbation. Both of these 
arguments established that the perturbation by a long 
range interaction is relevant for the short range spin glass 
if cr < cr,, = d. It is then recalled that the droplet theory 
predicted that the critical value cr^ depends on the stiff- 
ness exponent for the droplets as: (Jc = d—6. The combi- 
nation of these two exponent relations suggests that the 
droplet theory itself predicts that the stiffness exponent 
for droplets is zero: 9 — 0. 

While our analytical considerations are self contained 
and do not require input or justification from numeri- 
cal simulations, it is nevertheless reassuring that numer- 
ical results are in agreement with our conclusions. Katz- 
graber and Young [20|, for example, numerically studied 
the spin glass problem in d = 1 with a long range in- 
teraction with exponent cr. Working at cr = 0.75, they 
determined the droplet stiffness exponent 9 with great 
accuracy. The droplet theory of [1| predicted a droplet 
stiffness exponent 9 = d — a = 0.25 for this case. Katz- 
graber and Young, however, found 9 = —0.005, clearly 
inconsistent with the droplet prediction, but consistent 
with our analytical result of = 0, as well as with the 
TNT picture. It is further noted that Ref. [lOl also stud- 
ied the domain wall exponent and found it consistent 
with the TNT picture. The very recent work by Leuzzi 
et al is also indicative of 6* = 0.[13[ 

We also mention that P(A) was explicitly analyzed nu- 
merically in a one dimensional version of the above long 
range interacting model in a related recent paper. [21 1 
While P{X) exhibited a strong suppression approaching 
A = for cr < 1, the data did not compel the conclu- 
sion of P(A = 0) = 0. However, we believe that further 
work is in order here, as measuring a small exponent in 
the presence of large finite size effects requires a careful 
analysis. Our ongoing work on this issue will be reported 
elsewhere. [2^ 

The stiffness exponent being zero also has support 
from papers directly simulating short range spin glasses. 
Kisker ct al. studied the coarsening dynamics of the 
3d Edwards- Anderson Spin Glass and concluded that 
the scaling analysis of the autocorrelation function and 
correlation length in particular were suggestive of the 
stiffness exponent 9 = O.jl^ Their data were best 
fit with a droplet energy cost which depended loga- 
rithmically on the linear size: Edropiet{L) ^ ln[L). 
The vanishing of the droplet stiffness exponent 9 was 
also concluded from equilibrium studies of system size 
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excitations, [i, [13, [Til, [H, [H, [H | 

One could propose to resolve the above conundrum b 
arguing that the original exponent relation derived in 
was only a lowest order calculation and it needs to be 
reanalyzed in light of the importance of the avalanches. 
However, this only proves our point once again. The 
calculation in [^ was carried out assuming that the sys- 
tem was stable against perturbations. In stable phases, 
such lowest order perturbative calculations should be ap- 
propriate. Therefore, if it is argued that the avalanches 
rearrange the state of the short range spin glass to such 
a degree that a higher order, or even a re-summcd in- 
finite order perturbation calculation is necessitated be- 
cause e.g. P(A) develops a soft gap, that is essentially 
a recognition that the long range perturbations were in- 
deed relevant, proving our point once more. 

Before closing we mention that no explicit analysis was 
given above for the domain wall stiffness. Deriving the 
analogous exponent relation and the consequences of the 
pseudogap formation for domain walls is the target of 
ongoing work. 

In conclusion, wc have analyzed the stiffness exponent 
6 of the short range Ising Spin Glass by applying the 
probe of an additive perturbative long range interaction. 
We found that the long range interactions induce large 
avalanches in the system because of the large number of 
low-stability spins in the short range Ising Spin Glass. 
These avalanches roll until the system develops a pseudo 
gap at low energies, characterized by a power law distri- 
bution of the stabilities P(A) ~ A", where a = {d/a) — 1 
for d/2 < a < d and a = 1 for cr < d/2. This large 
scale reorganization makes the long range interaction a 
relevant perturbation for all a < ac ~ d exponent values. 
Finally, comparing this critical exponent with the predic- 
tion of the droplet theory oi ac = d — 6, one concludes 
that combining the logic of the droplet theory with the 
notions of avalanches and pscduogaps points to the stiff- 
ness exponent being 9 = in the short range Ising Spin 
Glass. While the vanishing of the droplet stiffness has 
consequences for the droplet theory, it is consistent with 
the RSB and TNT theories, both of which are consistent 
with such a weak stiffness. 
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